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ABSTRACT 

We present results of the evolution of AGB stars of 3 Mq and 5 Mq with solar metallic- 
ity calculated with the Eggleton stellar evolution code (stars), which has a fully im- 
plicit and simultaneous method for solving for the stellar structure, convective mixing 
and nuclear burning. We introduce the concept of a viscous mesh in order to improve 
the numerical stability of the calculations. For the 5Mq star, we evolve through 25 
thermal pulses and their associated third dredge-up events. We obtain a maximum 
helium luminosity of 1.7 x 10^ Lq and significantly deep dredge-up after the second 
pulse. Strong hot-bottom burning is observed after the 5"^ pulse. The 3 Mq model is 
evolved through 20 thermal pulse events and we find third dredge-up after the 7th 
pulse. During the 14th pulse sufficient carbon has been brought to the surface to pro- 
duce a carbon star. We find that dredge-up and the transformation into a carbon star 
occur at significantly smaller core masses (0.584 and 0.608 M©, respectively) than 
in previous calculations for 3M0. 
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1 INTRODUCTION 

. The observed enhancements of carbon and s-process el- 
' ements in asymptotic giant branch (AGB) stars show 
that material, processed by helium burning during ther- 
mal pulses, is being dredged up to the surface. Although 
qualitatively well understood, the extent and efficiency of 
this third dredge-up (TDUP) process has been very hard to 
model quantitatively, and remains one of the major uncer- 
tainties in AGB evolution. The amount of TDUP obtained 
for a star of a certain mass and composition differs greatly 
between one model calculation and another. It depends on 
many factor s, such as the numerical treatment of convective 
boundaries iLatta nzio 198S H, whether allowance for extra 
mixing IS made (|H o11owc11 fc Ibenlll98^. the mixing-length 
parameter iBoothrovd fc Sackmannll^88) and the numeri- 
cal re solution of the mass grid and time step dStraniero et alJ 
Il997t) . Recently, very efflcient dredge-up for stars with small 
core masses has been found by various authors (Herwig et 
al. 1997, Mowlavi 1999, Herwig 2000) if a diffusive form of 
convective overshooting is taken into account. 

In the calculation of the evolution of stars of 1 — 
8M0 through the thermally pulsing asymptotic giant 
branch (TP-AGB) phase of their evolution most evo- 
lutionary codes use only a partially self-consistent ap- 
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proach. The treatment of structure, mixing and nucle- 
osynthesis is not simultaneous, with mixing often be- 
ing solved in a separate iteratio n step from th e othe r 
two. Examples include the co d es of Straniero et al.l lll997tl . 
Wagenhuber fc Groeneweeenl il99g ). iHerwid ll2000h and 
Karakas^attami^m^Poy (12002*) . As the phenomenon of 
third dredge-up depends critically on the treatment of con- 
vection within a stellar structure calculation it is desirable 
to combine the calculation of nucleosynthesis, mixing and 
structure into a single, s imultaneous s tep. T he first attempt 
to do this was made bv lPols fc Tbu3 teOOlll . In their calcu- 
lations a 5 M0 star was evolved without mass loss through 
the first 6 thermal pulses and deep dredge-up occured. 

The occurence of deep dredge-up is important for the 
formation of carbon stars. These are defined as stars that 
are M-type and have surface carbon-to-oxygen abundance 
ratios exceeding unity. Observational evidence suggests that 
these stars are of low mass, most likely between 1 and 3 M0 
(Iben 1981). However, there has been considerable difficulty 
in producing detailed theoretical models of carbon stars with 
low enough masses and luminosities. Of the early work on 
the subject, iBoothrovd fc SackmannI Jl988D were able to 
produce two carbon star models from initial masses of 1.2 
and 2.0 Mq und er metal-poor conditions {Z — 0.001), while 
iLattanziol il989() produced a model initially of M = 1.5 Mq 
with Z = 0.02. More rece ntly low-mass carbon star mod- 
els have been produced bv lStraniero et al] (Il997f) and, with 
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the aid of convective overshooting. iHerwid i200Ci) . However, 
the core masses found in these models are still too large to 
explain the carbon star luminosity function (Izzard & Tout 
2004). 

While the approach of Pols & Tout enabled them to 
make significant progress with a fully implicit and simulta- 
neous treatment of TP-AGB calculations, problems with nu- 
merical stability prevented them from evolving their model 
further than a few thermal pulses. In this work we have 
developed methods that overcome some of the numerical in- 
stability that occurs in a fully implicit treatment of the evo- 
lution. These methods are described in section 2. Results of 
the calculations for a 5Mq star are presented in section 3 
and compared to other available models in section 4. In sec- 
tion 5 we present results of the calculations for a 3 Mq and 
detailed comparisons are made to other models in section 6. 



present and previous timestep, At is the timestep, Sm^._j_i 
are the masses contained in the zones above and below mesh- 
point k, Rx,k is the net consumption rate of X by nuclear 
reactions, and (t^,^ j. are the difi^usion coefficients correspond- 
ing to these zones. The linear diffusion coefficient D is re- 
lated to a by cr = {iirr p) D, where r and p are radius and 
density. 

Convective mixing is is treated in the f ramework of 
the mixing- length theory jBohm- VitensdllQSd) . which gives 
D — DyiuT ~ ^vl with V and I the mean velocity and mean 
free path of convective eddies. In order to maintain numer- 
ical stability, iPols fc Touli 0200 li') found it useful to express 
D as 



D = Dmlt 



with 



1~{1-I3)W 



(3) 



2 THE STELLAR EVOLUTION CODE 

The models were calculated using the STARS ID st ellar evo- 
lution code originally developed bv lEggletonI lll97J) and up- 
dated by many contributors (e.g. Pols et. al. 1995). This code 
solves the structural and composition equations, including 
the effects of both mixing and burning, implicitly and si- 
multaneously. We compute the structure of a star over a 
given number of mesh points which are placed in the most 
appropriate regions determined by equal intervals of a mesh 
spacing function Q{P, T, r, m) where P, T, r and m are pres- 
sure, temperature, radius and mass respectively (lEggletonl 
I197J) . This function is designed to push more mesh points 
into regions where they are needed such as burning shells 
and ionization zones. The mesh is adaptive so that as the 
star evolves, the mesh points move around to stay with phys- 
ically important regions. We apply a mesh-spacing function 
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where ci . . . cg are appropriately chosen constants, M is the 
total mass and Ph and Phc are the values of P at the posi- 
tion of the H- and He-exhausted core mass, respectively. The 
last two terms are included only during the TP-AGB phase. 
They concentrate a large part of the mesh, about three- 
fifths, in the intershell region (between O.IPh and 3Phc) 
while at the same time ensuring that during a TP, when 
T and P undergo rapid changes, the mesh does not move 
around too much in the intershell region and so that numer- 
ical diffusion is suppressed. The mesh-spacing equation is 
also solved implicitly and simultaneously with the structure 
and composition equations. 

The code includes mixing and nuclear burning using a 
diffusion equation approximated by an implicit second-order 
difference equation of the form 
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The value of /3 was treated as a free parameter and allowed 
to take on different values in the H-rich envelope and in the 
H-free intershell region and core. 

All the quantities including o are defined only at each 
meshpoint k so it is not a priori obvious how o^^^i should 

be calculated. In our standard calculations we take (Ji.^i = 

1 
2 



where Xk and X° are the abundances at meshpoint k at the 



(cTfc-hafc+i), following |EggletonHl97d) . |Pols fc Toutl^l2001^ 

investigated the effect of taking a geometric mean = 
ycTfe ■ Ok+i rather than an arithmetic mean. These schemes 
effectively differ only in zones where the Schwarzschild 
boundary is located between meshpoints k and fc -I- 1. The 
arithmetic mean in effect allows the radiative meshpoint ad- 
jacent to a convective boundary to be mixed while the geo- 
metric mean does not. The choice is usually inconsequential 
but becomes important in the presence of a composition 
discontinuity that leads to a discontinuity in Vrad through 
the opacity dependence. Such discontinuities arises during 
core-helium burning at the edge of the convective core and 
at the bottom of the convective envelope during TDUP 
(Paczyhski 1977). The standard arithmetic mean ensures 
that Vrad " Vad always approaches zero from the convec- 
tive side of a boundary, even in the case of a discontinu- 
ity, so a convective boundary always corresponds to a stable 
Schwarzschild boundary (Vrad = Vad) - However when there 
is a discontinuity in Vrad at the boundary even the slight- 
est extra mixing makes it unstable and physically we expect 
mixing to extend until Vrad ^ Vad when material can mix 
across the boundary. Pols and Tout demonstrated that use 
of the geometric mean for o-j,_^i suppresses this physical be- 
haviour and prevents third dredge-up (see also the discussion 
in Mowlavi 1999 and Herwig 2000) so we always employ the 
arithmetic mean. 

The TP-AGB is a difficult phase of evolution to calcu- 
late. As the peak of a thermal pulse is reached the code must 
use very small timesteps, down to the order of a few minutes. 
This can cause problems with the stability of the code which 
we outline below. After a pulse has occured it is necessary 
to limit the timestep to the order of weeks in order to prop- 
erly resolve the process of TDUP. After this the code needs 
to be able to increase its timestep significantly to deal with 
the long quiescent interpulse period before the next thermal 
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pulse. Sections 2.1 and 2.2 describe in detail the numerical 
measures we have added with these constraints in mind. 

2.1 The Viscous Mesh 

At short timesteps, such as those required to resolve the peak 
helium luminosity of a thermal pulse, there is a numerical 
instability associated with the luminosity equation 

Lk+i - Lk = {m'Ei)f^^i + {m'E2)k[mk] 

-{m'E2)k+i[-mk+i] (4) 

where Lk is the luminosity at mesh point number k. The 
change in mass with mesh point is denoted m' = dm/dk. 
Square brackets signify a term only included when it is pos- 
itive. The first term containing Ei is related to the usual en- 
ergy generation terms but evaluated at constant mesh point 
and the terms containing E2 are an upstream approxima- 
tion for the advection term owing to the adaptation of the 
mesh. When the timestep becomes small these terms become 
large and this leads to numerical instability when two large 
numbers are subtracted to give a small result. 

To deal with this the mesh can be fixed in mass and 
so eliminate the last two terms of the above equation. How- 
ever, it is not desirable to fix the entire mesh because the 
outer regions do not cause the same problems as the inner 
ones and the mesh need only to be fixed at small timesteps. 
This causes a second problem - how do we return from the 
fixed to the fully adaptive mesh? The mesh cannot just be 
instantaneously returned to a fully adaptive state because, 
during the fixed phase, mesh points drift away from where 
the mesh spacing function would place them. Rather, the 
mesh must be returned to a fully adaptive state gradually 
and so we have developed the idea of a viscous mesh. 

Our viscous mesh combines fixed and adaptive be- 
haviours. A weighting coefficient 7, a function of mesh point 
number and the timestep size, is used to determine the na- 
ture of the mesh. In problematic areas the mesh is gradually 
fixed by increasing the value of 7. If 7 = 1 a mesh point is 
fully fixed in mass while if 7 = the point is fully adaptive. 
We implement this by solving 




alongside the equations of stellar structure and composition. 
We choose the weighting coefficient 7 to be a function of the 
timestep with 7 becoming unity for timesteps of 10~^ yrs 
and below for the central mesh points. 



2.2 Timestep control 

The STARS code attempts to choose the most appropriate 
timestep size based on the changes to the variables of the 
previous model required to produce the current model. A 
sum of the absolute values of these changes (excluding those 
made to the luminosity) is made over all variables and over 
all meshpoints producing a single numerical value 

d = EiEfe|a;i,fc|, 

where Xi^k are the values of the changes made to each vari- 
able i at a given meshpoint k. Usually the value of d is 
dominated by the temperature and the degeneracy. This is 



then compared to a preset optimum value dopt. If dopt/d 
is greater than one the timestep is increased by this frac- 
tion or 1.2, whichever is smaller. If dopt/d is less than one 
the timestep is reduced by this fraction or 0.8, whichever is 
larger. In this way the most appropriate timestep is chosen 
for the next model. 

The value chosen for dopt depends on the number of 
mesh points in the model. For a 999 mesh point model, a 
value of 5 is typically chosen. However, if it is desirable to 
have smaller timesteps then lower values are used. In order 
to deal with the range of timescales involved in TP-AGB 
evolution we find it necessary to change the value of dopt at 
certain stages of the evolution. 

As the peak of a thermal pulse is approached the 
timestep drops to about 2 x 10~^ years. To avoid numerical 
problems inherent at lower timesteps we need to boost dopt. 
We typically do this by adding 0.5 when the timestep drops 
below 2.5 X 10~^ years. When the helium luminosity begins 
to decline after the peak of the thermal pulse the intershell 
convection zone (ICZ) shuts down and the convective en- 
velope of the star moves inward. During this phase dopt is 
slowly returned to 5. As the envelope reaches down into the 
region where the ICZ was active it begins to dredge-up the 
products of helium burning. At this point it is necessary to 
limit the size of the changes being made to the variables to 
avoid numerical instability. 

We have written an algorithm that records the position 
of the boundary of the ICZ at its maximum outward extent. 
We then compare this to the location of the boundary of 
the convective envelope. When the envelope reaches into the 
region where the ICZ was active the value of dopt is set to 1. 
This limits the timestep to the order of a few days. When the 
helium luminosity reaches as low as 3 x 10"^ L0 and TDUP 
is over dopt is restored to 5. 



3 CALCULATIONS FOR A 5 STAR 

We apply our code to a 5 M0model with Z = 0.02, evolved 
from the pre-main sequence with 199 mesh points and so- 
lar composition. Before the onset of the TP-AGB phase we 
increase the number of mesh points to 999 to provide suf- 
ficient spatial resolution in the intershell region. Through- 
out the evolution we do not consider mass loss nor do we 
include convective o vershooting. We use the formalism of 
iBohm-Vitensd il958h for convection with the ratio between 
the mixing length and the pressure scale height a = 1.925. 
We choose our /3 = 1 in the envelope of the star, and 
/3 = 5 X 10~^ in regions where the ma ss fraction of hydro - 
gen is less than lO"''. As determined bv lPols fc Tout! (1200 ill 
there is negligble quantitive difference if we increase these 
values. 

In the course of these evolutionary calculations prob- 
lems were encountered with the composition equations. 
These were found to be related to the spontaneous creation 
of convective zones above and below the ICZ. An exam- 
ple of the lower zone is shown in Figure Q The zones be- 
low the ICZ occur in regions where the temperature gra- 
dient is positive and hence are not true convective regions. 
At most these zones are no more than a couple of mesh 
points wide and their lack of extent confirms that they are 
not failed degenerate pulses as described in the work of 
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Model Number 



Figure 1. Plot showing a spurious convection zone below the 
helium burning shell. The boundaries of convection zones are 
marked with solid lines. The hydrogen and helium burning shells 
are marked with short and long dashed lines respectively. The 
zone occurs in a region where the temperature gradient is posi- 
tive and so does not represent a true convective region. The upper 
zones cannot be seen in this figure as they lie between the top of 
the ICZ and the hydrogen burning shell. 



iFrost. Lattanzio fc WoodI ijlQQST) . These zones begin to oc- 
cur after the 7"^ thermal pulse but do not cause serious 
problems with model convergence until the 12"^ pulse. 

In order to eliminate these spurious zones we apply the 
following strategy. At timesteps below 4 x lO"** years we 
reduce the value of f3 to 10~^ for the 25 mesh points im- 
mediately above the top of the ICZ. We do this because we 
expect this zone to be radiative and so no mixing should oc- 
cur there. The spurious zones do not then appear. Below the 
helium burning shell we prevent the use of both convective 
mixing and the transport of energy by convection. These 
measures have been tested both on the 6 pulses where no 
spurious zones were found and on pulses 7 to 12 where they 
occurred. The evolution is the same both with and with- 
out them. Thus we are confident that they do not affect the 
evolutionary calculations. 

With the outlined modifications, a total of 25 thermal 
pulses and their associated TDUP episodes were calculated. 
A plot of the evolution of the helium luminosity with time 
is shown in Figure |5| and the evolution of the hydrogen- 
and helium-exhausted core masses with time is shown in 
Figure 0] These core masses are defined as the mass coor- 
dinates where the mass fractions of H and He are less than 
0.3. Thermal pulses begin when the helium core mass is 
0.838 M0. We find the peak helium luminosity rapidly ap- 
proaches 1.7 X 10^ L0. Third dredge-up occurs on the very 
first pulse and significant hot-bottom burning is seen after 
the 4*'' pulse. 

Owing to the efficiency of third dredge-up we find that 
the intershell region narrows initially, as shown in Figure |21 
When the phase of most efficient dredge-up is over the inter- 
shell continues to narrow because the helium burning shell 
moves outward at a greater rate. This is to be expected be- 
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Figure 2. Plot of the evolution of the helium luminosity versus 
time since the peak of the first thermal pulse in years, for the 
5M0 model. 
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Figure 4. The evolution of the hydrogen-exhausted (solid line) 
and helium-exhausted (dashed line) core masses with time since 
the peak of the first thermal pulse, for the 5Mq model. 



cause the helium luminosity in the interpulse phase is greater 
in the later pulses. 



3.1 The Effects of TDUP 

The efficiency of TDUP is usually expressed in terms of a pa- 
rameter A which is the ratio of the amount of mass dredged 
up (AA'/dup) to the amount by which the hydrogen ex- 
hausted core grows in the preceeding interpulse (AMc), i.e. 
A = AA/dup/AMc. The value of A for each pulse is shown in 
Table We find that A quickly reache s a value of unity b y 
pulse 5, in agreement with the work of IPols fc Tout! i200ll) . 
We then find that it slightly exceeds unity for a number of 
pulses, before dropping back down to just below it. It should 
be noted that the A values obtained h ere for the fi r st few 
pulses are greater than those found bv lPols fc Tout! (|20^) 
refiecting improvements in the accuracy of the solution of 
the equations in the core. 

The evolution of the surface abundances of ^^C, ^''N 
and ^''O in our model is shown in Figure]^ Initially, the car- 
bon abundance increases as TDUP draws carbon produced 
by helium burning to the surface. As dredge-up becomes 
deeper the envelope reaches down into hotter and hotter 
layers. Eventually it reaches a depth where the temperature 
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Figure 3. Plot of the structure of the intershcU region versus model number. The convective zone boundaries are displayed as solid 
lines. The hydrogen and helium burning shells are shown by short and long dashed lines respectively. Initially, efficient third dredge-up 
causes the intershell region to diminish. Over the latter pulses, the helium burning shell moves outward at a greater rate. We use model 
number rather than time in order to resolve the structure of the pulses. The ICZ typically exists for 20 years while the interpulse period 
is about 13,000 years. 
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Table 1. Details of our 5MQmodel. The data are: TP - the thermal pulse number, A/jj - the hydrogen free core mass, Tip - the interpulse 
period, LJJ^" - the peak luminosity from helium burning, AAfjj - the hydrogen free core mass growth during the interpulse, AA^dup - 
the amount of material dredged up, A the dredge- up efficiency parameter and C/O - the surface carbon-to-oxygen ratio by number. 
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Figure 5. Surface abundance evolution with time for '^^C,'^*N 
and for the 5Mq model. While TDUP draws fresh carbon 
to the surface it is destroyed by hot-bottom burning during the 
interpulse. 



is high enough for the CNO cycle to begin. This is hot- 
bottom burning (HBB). The efTect of HBB is to deplete the 
surface of carbon and oxygen whilst enhancing the levels of 
nitrogen. By the 5*'^ thermal pulse HBB is strong enough to 
remove more carbon than is dredged up. Hence our models 
show a drop in the surface carbon abundance. 

Once the phase of deepest dredge-up is passed we reach 
a phase where carbon enhancement by TDUP is slightly 
greater than carbon depletion by hot-bottom burning. We 
find a small but persistent increase in the surface carbon 
abundance. The persistent rise in the nitrogen surface abun- 



dace is due to hot-bottom burning and by pulse 11 it is the 
most abundant of the CNO elements. The oxygen level is 
only slowly depeleted because the temperature of the base 
of the convective envelope is not high enough for the NO 
part of the cycle to operate efficiently. 



4 MODEL COMPARISONS 

Various groups have produced models of a 5Mq star with 
Z=0.02 and no mas s loss. We compare the ir results to ours. 

The model of IStraniero et alJ (|200(J) . calculated us- 
ing the Frascati Raph son-Newton Evolutionary Code 
JChieffi fc Stranierd igS?). produces a peak helium luminos- 
ity of 1.3 X 10® Lq . This is lower than our value by an order of 
magnitude. Their model also approaches the peak luminos- 
ity much more slowly, taking around 20 thermal pulses to get 
there. We reach peak luminosity after just 10 pulses. As we 
discuss in Section 6, we suspect that the higher peak lumi- 
nosity is related to the deep dredge-up we find, although we 
cannot verify this suspicion in this case because the dredge- 
up efficiency of their model is not available in the literature. 
Interestingly Straniero et al. do not find hot-bottom burn- 
ing in their model. They give the temperature at the base 
of the convective envelope as 4.6 x 10^ K which is not hot 
enough for the CNO cycle to be effective. Our model has 
a temperature of 7.9 x 10^ K at the base of the envelope 
owing to penetration down to the deeper, hotter regions of 
the star. 

In comparison t o the model of 

iKarakas. Lattanzio and Polj l|20o3), produced with the 
Monash version of the Mount Stromlo Stellar Structure 
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Figure 6. Plot of the evolution of the hydrogen-exhausted (solid 
line) and helium-exhausted (dashed line) core masses with time 
since the peak of the first thermal pulse, for the 3 Mq model. 
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Figure 7. Plot of the evolution of the luminosity from helium 
burning with time since the peak of the first thermal pulse, for 
the 3 Mq model. 



Program (MSSSP), we find a significantly smaller core mass 
at the first thermal pulse. Karakas et al. give a core mass 
of 0.862 M© whereas we calculate the value to be 0.838 Mq. 
Comparing the efficiency of third dredge-up between the 
two models we find ours to be more efficient. They find 
Amax of 0.955 for their equivalent model, compared to our 
value of 1.05. 



5 CALCULATIONS FOR A 3 Mq STAR 

The effects of hot-bottom burning in the 5 Mq model pre- 
vent the build-up of carbon on the surface of the star. 
Stars of lower mass do not have envelope temperatures high 
enough for HBB to occur and so there is a possibility that 
a carbon star could be formed. 

We apply our code to a 3 Mq star evolved from the 
pre-main sequence with 199 mesh points. At core helium 
burning the number of mesh points used is increased to 499 
in order to facilitate the transition to a model with 999 mesh 
points when the TP-AGB is reached. We choose a = 1.925 
and ignore mass loss and convective overshooting, as for the 
5Mq case. 

We evolve a total of 20 thermal pulse events. TDUP 
occurs after the 7th pulse when the hydrogen-free core mass 
is 0.584 Mq. A plot of the evolution of the hydrogen- and 
helium-free core masses is shown in Figure|S| The C/0 ratio 
exceeds unity during TDUP following the 14th pulse at a 
core mass of 0.608 Mq. The pulses reach a maximum helium 
luminosity of Lhc = 5.9 x 10* Lq. A plot of the evolution 
of the helium luminosity with time is shown in Figure |7| 

The characteristics of the model are listed in Table |21 
and the effects of TDUP on the surface abundances by mass 
of the CNO elements can be seen in Figure|H| Initial dredge- 
up episodes are quite weak but the efficiency rises rapidly 
at each of the first few dredge-up events. Large amounts of 
carbon and smaller quantities of oxygen are carried to the 
surface. By pulse 20 dredge-up is deep with A = 0.901. 



6 DETAILED MODEL COMPARISONS 

Many groups have produced models with similar initial con- 
ditions to ours. Here we compare the characteristics of these 
models. 



0.016 
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Figure 8. Plot of the evolution of the CNO element abundances 
by mass with time since the peak of the first thermal pulse, for 
the 3 Mq model. 

The unique aspect of the STARS code is that it is fully 
self-consistent - there is no splitting into separate steps of 
the solution for the equations governing structure, nuclear 
burning and mixing of elements. Most codes solve for struc- 
ture and burning together and then apply the mixing in 
a separate step based on this solution. This approach is 
not guaranteed to yield a self-consistent solution. Exam- 
pl es of models produced by t his method include the work 
oflStraniero et all (119971) and iKarakas. Lattanzio and Pold 
i2002h . Both have produced 3 Mq models at solar metallic- 
ity. 

The models of IStraniero et al.l ^1997^ were made us- 
ing the Frascati Raphson Newton Evolutionary Code 
iChieffi fc Stranierdil989.,) . They have evolved models both 
with and without mass loss. Details of the model without 
mass loss are given in Table |3 While the model has not 
been evolved to the point where it will become a carbon 
star, it is very likely to do so within the next few thermal 
pulses. 

Straniero et. al.'s model starts of with a slightly more 
massive core, 0.007 Mq heavier than ours. They find that 
dredge-up does not start until the core mass has reached 
0.611 Mq. In contrast we find that TDUP begins at a core 
mass of 0.584 Mq . Our model also shows a much deeper level 
of dredge-up. The peak value of A in Straniero et al's model 
is 0.461 compared with a value of 0.903 in our model. We 
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Table 2. Details of our 3 Mq model. The data are: TP - the thermal pulse number, A/jj - the hydrogen free core mass, rip - the interpulse 
period, ^^^^ - the peak luminosity from helium burning, AAfjj - the hydrogen free core mass growth during the interpulse, AAfjjup - 
the amount of material dredged up, A the dredge-up efficiency parameter and C/O - the surface carbon-to-oxygen ratio by number. 
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Table 3. Details of a 3M0 star evolved without mass-loss, taken from Table 1 of IStraniero et all il997l) . The data are: TP - the 
thermal pulse number, Mjj - the hydrogen free core mass, rjp - the interpulse period, AAfpj - the hydrogen free core mass growth during 
the interpulse, log(Lg^^''/ Lq - the peak luminosity from helium burning, A the dredge-up efficiency parameter and C/O - the surface 
carbon-to-oxygen ratio. 



also note a difference in the behaviour of the peak helium 
luminosity. The pulses of our model are more luminous and 
the luminosity increases more rapidly from pulse to pulse. 

The phenomena of deep dredge-up and high peak lumi- 
nosities are related, as was also found by Straniero et al. On 
the one hand, stronger shell flashes cause greater expansion 
during the power-down phase and thereby allow the con- 



vective envelope to reach deeper layers. Conversely, deeper 
dredge-up leads to a longer interpulse period and greater val- 
ues of AMh as found in our model (Table 2). Hence a larger 
reservoir of helium is built up which may cause a larger 
peak helium luminosity. Another factor that contributes to 
the more powerful shell flashes is that the longer interpulse 
period gives the inactive He-burning shell more time to cool 
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and be compressed. Therefore the next shell flash occurs at 
higher density and more work has to be done to expand the 
intershell region, i.e. a higher peak luminosity is required. 
The upshot is that deep dredge-up and powerful shell flashes 
mutually reinforce each other. 

In comparision w it h the models of 
iKarakas. Lattanzio and Polj i2002h . calculated with 
the Mount Stromlo Stellar Structure Program, we see 
similar trends. Their model has a core mass of 0.595 Mq 
at the first thermal pulse and TDUP begins when the core 
mass reaches 0.635 M©. These core masses are significantly 
heavier than those obtained in our model. They find a 
maximum dredge-up efficiency of 0.790 which is also lower 
than the value we obtain. 



6.1 Models with Convective Overshooting 

It is generally agreed that models produced using only con- 
vective mixing do not produce the correct abundance profiles 
for s-process nucleosynthesis to occur. They do not produce 
a pocket of ^^C in the intershell region which is required 
to provide a neutron source via the reaction ^^C(Qf,n)^'^0. 
An additional mixing mechanism is required and no physical 
origin has yet been established. One of the mechanisms that 
have been proposed is convective overshooting, where blobs 
of material penetrate into convectively stable regions leading 
to partial mixing of protons and carbon-rich material. 

TP-AGB mo dels wi t h con vective overshooting have 
been published bv lHerwi3 i2000h . He has calculated the evo- 
lution of both 3 and 4 Mq stars in which he allows convective 
elements to penetrate the convectively stable regime with an 
exponentially decaying velocity. This overshooting is applied 
at all convective boundaries. These models are evolved with 
a Reimer's mass-loss law 

r r> 

M = -4x lO-i^:^^M0yr-' 

where L, R and A/ are the luminosity, radius and mass of 
the star in solar units. 

Selected details of Herwig's model, taken from Table 1 
of lHerwia(.2000i) . are listed in Table |1| Note that the values 
of A obtained are similar to our calculations even though 
we do not include overshooting in our models though our 
dredge-up efficiencies do not exceed unity. 

This model naturally starts off with a higher core mass 
than our model because the evolution prior to the AGB 
with convective overshooting produces larger cores. Accord- 
ing to Herwig the large dredge-up efficiency in his model 
is the result of overshooting below the pulse-driven inter- 
shell convection zone. This leads to a higher temperature 
at the bottom of this zone and a larger peak luminosity 
than without overshooting. Stronger shell flashes give rise 
to deeper dredge-up, as discussed above. In contrast we find 
similarly high peak luminosities and similarly deep dredge- 
up even without overshooting. However, it should be noted 
that Herwig's models produce a ^'^C pocket (albeit one that 
is too narrow - see Denissenkov & Tout 2003) whereas ours 
require another mechanism. 
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Table 4. Selected details of a 3 Mq star evolved with convectiv e 
overshooting and mass- loss taken from Table 1 of lHerwid i200(t . 
The data are: TP - the thermal pulse number, LHe - the peak 
luminosity from helium burning, A the dredge-up efficiency pa- 
rameter, Mjj - the hydrogen free core mass and C/O - the surface 
carbon-to-oxygen ratio. 



7 CONCLUSION 

We have made a fully implicit calculation, with simultane- 
ous solution of the structure, mixing and nuclear burning, of 
the evolution of both 5 Mq and 3 Mq stars. We have evolved 
the 5 M0 star through 25 thermal pulses and their associ- 
ated third dredge-up events. We find much stronger peak 
helium luminosities than are obtained by non- simultaneous 
codes and more efficient third dredge-up. Hot-bottom burn- 
ing is found to occur during the interpulse. This prevents 
the build-up of carbon at the surface and leads to nitrogen 
being the most abundant CNO element. 

We have followed the evolution of a 3 Mq star through 
20 thermal pulses. We find deep third dredge-up to occur 
after the 7'*' thermal pulse and the star becomes a car- 
bon star after the 14**^ thermal pulse. We find that these 
events occur at much lower core masses than those found 
by calculations with non- simultaneous codes. We also find 
much more efficient third dredge-up. The presence of effi- 
cient dredge-up and low core mass may help to explain the 
carbon star luminosity functi o n. Wo rking with the results of 
iKarakas. Lattanzio and Polj ^00^, Izzard & Tout (2004) 
demonstrated that, in order to fit the observed LMC and 
SMC carbon star luminousity functions, the minimum core 
mass for dredge-up to occur must be 0.07 Mq lower than 
those calculated with the MSSSP. While our models are at 
higher metallicity, we find dredge-up to occur when the core 
mass is 0.584 Mq or 0.05 Mq lower than the corresponding 
model of iKarakas. Lattanzio and Polj (|20o3). 

We conclude that our calculations yield results that dif- 
fer in very important ways from calculations made with non- 
simultaneous codes. We find third dredge-up to be more ef- 
ficient and that it occurs at lower core masses. This may be 
due to the simultaneous solution of the equations although 
there are other minor differences between the numerical and 
physical details of different evolution codes. It should be 
noted that we have not conisdered the effects of mass loss in 
this work. Third dredge-up depends on the envelope mass 
(see e.g. Straniero et al. 2003) and the inclusion of mass loss 
would be expect to reduce the efficiency of third dredge-up. 
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